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3 Abstract. The spatial and spectral characteristics of the turbulent plasma 

4 density, electric fields, and ion drift in ionospheric E-region are studied us- 

5 ing a new set of nonlinear plasma fluid equations. The fluid model combines 

6 both Farley-Buneman (Type-I) and Gradient-Drift (Type-II) plasma insta- 

7 bilities in the equatorial electrojet. In our unified model of the plasma in- 
a stabilities we include the ion viscosity in the ion momentum equation and 

9 electron inertia in the electron momentum equation. These two terms play 

10 an important role in stabilizing the growing modes in the linear regime and 

11 in driving the Farley-Buneman instability into the saturation state. The sim- 

12 ulation results show good agreements with a number of features of rocket 

13 and radar observations, such as: (1) saturation of plasma density perturba- 
« tions depends on the solar condition and reaches 7-15% relative to the back- 
is ground, (2) fluctuation of the horizontal secondary electric field reaches 8- 

i6 15 (mV/m), (3) stabilization of the phase velocity of the perturbed density 

I? wave around the value of the ion-acoustic speed inside the electrojet, (5) "up- 

18 down" asymmetry in the vertical fluxes of the plasma density, (6) "east-west" 

19 asymmetry of the plasma zonal drifts, and (7) generation of small-scale of 

20 the order of meter scale lengths irregularities embedded in large-scale struc- 

21 tures. Spectral analysis of the density fluctuations reveals the energy cascade 

22 due to the nonlinear coupling between structures of different scales. The break- 

23 up of the large-scale structures into small-scale structures explains the dis- 

24 appearance of Type-II echoes in the presence of Type-I instabilities. 
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1. Introduction 

25 The Earth's ionosphere is rich with plasma irregularities of various scale-lengths ex- 

26 tending from a few centimeters to hundreds of kilometers [Farley and Balsley 1973, Fejer 

27 1980]. At times during which the ionosphere is disturbed, mixed plasma states occur 
2a with large-scale coherent structures such as plumes, bubbles and vortices along with high 
29 levels of meter-scale incoherent plasma turbulence. These dynamics are important from 
so both fundamental and applied perspectives. The combination of small-scale turbulence 
si with large-scale coherent structures is at the forefront of basic plasma turbulence theory, 

32 while the spatial and spectral characteristics of the evolution of the electron density in 

33 the E-region is important for terrestrial navigation and communication systems [Wernik 

34 2003] which resolution is often limited by the turbulence in ionospheric plasma density. 

35 In the equatorial E-region of Earth's ionosphere, the Farley-Buneman instability is 

36 driven by an intense equatorial eastward electrojet current [Richmond 1973, Rishbeth 

37 1997]. The Farley-Buneman instability is identified in the backscattering radar data as 

38 Type-I turbulent echoes, e.g. as observed by the coherent backscattering radar at Jica- 

39 marca [Bowles et al. 1963, Fejer 1975, Hysell 2007] and theoretically studied by Farley 
« [1963] and Buneman [1963]. 

41 The electrons are magnetized [y en <C oj ce ) and experience -UxB drifts in the E-region, 

42 while the ions are unmagnetized (u in 3> uj ci ) so they are tied to the neutral gas and ex- 

43 perience U drifts, where U is the speed of the neutral wind. The E-region, extending 

44 between altitudes 90-120 km, is found to be the most conductive region in the ionosphere 

45 owing to the large Cowling conductivity, Oc = crp(l + cli/vp), which depends on the ra- 
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« tio between the Hall, an, and the Pedersen, ap, conductivities [Baker and Martyn 1953]. 

47 This large Hall-to-Pedersen conductivity ratio, which reaches its maximum value of 15- 

48 20 between 104-105 km, allows the electron drift velocity, ve, to exceed the ion-acoustic 

49 speed C s = yj (T e + •y i T i )/m i by a factor that is large enough to drive the Farley-Buneman 

50 instability. 

51 In addition to the Farley-Buneman instability, density gradients can also drive another 

52 type of instability in the ionosphere. The gradient-drift instability is identified in the 

53 backscattering radar data as Type-II turbulent echoes, e.g. as observed by the coher- 

54 ent backscattering radar at Jicamarca [Bowles et al. 1963, Fejer 1975, Hysell 2007]. 

55 The density gradient drive instability mechanism was studied by Simon [1963] and Hoh 
se [1963]. From the bottom of the E-region to the peak of the F-region, the ionosphere 
s? is subject to strong drift wave type plasma turbulence owing to the upward increase of 
sa the electron density. The electron density between 90-110 km has an exponential profile 

59 n(z) = n 0 exp(z/L n ) perpendicular to the flowing current, y, and the geomagnetic field x, 

60 where z is the unit vector in the vertical upward direction and L" 1 = d z £nn is the electron 

61 density vertical gradient length scale. The gradient-drift instability can be excited by little 

62 or no threshold drift speed, which makes it more easy to excite than the Farley-Buneman 

63 instability. Farely and Balsley [1973] used the linear theory to get a relation between 

64 the threshold drift speed of the gradient-drift waves and their wavelength and have to 

65 be v E > k 2 Cl 1 + VenU ™/ UJ < xU ™ L Hence, we can have gradient-drift waves with a broad 

66 range of wavelengths; as long as kilometers and as small as tens of meters. However, it is 

67 not possible to find linearly growing gradient-drift waves of the order of meters because 

68 they require drift speeds much higher than what is measured and found in the equatorial 
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eg electrojet region. 

70 During the nighttime the recombination process happens very fast and we have regions 

?! of positive and negative gradient-drifts distributed randomly around the globe [Kudeki 

72 et al. 1982]. The gradient-drift instability is driven by the ExB drift instability due to 

73 gradients in electron density and dominates at small wavenumbers. 

74 Both Type-I and Type-II have been extensively studied in the past by many authors. 

75 In contrast, in this work, we introduce a unified model describing both the Type-I and 

76 Type-II modes in the ionospheric E-region. The dominant nonlinear coupling mechanism 

77 is the ExB convection of the fluctuating density described mathematically by a Pois- 

78 son bracket or a Jacobian nonlinearity [Horton and Hasegawa 1994, Horton and Ichikawa 

79 1997]. 

so This paper is organized as follows. In next section we review the rocket measurements 

si and radar observations and related theoretical work. In section(3) we present the physical 

82 model that describes the nonlinear dynamics of both Farley-Buneman and gradient-drift 

aa instabilities. The linear limit of this dynamical system and linear results are discussed in 

84 section(4). The nonlinear simulations and results are presented in section(5). We end the 

as paper with a summary of all results and the conclusions we can draw from these results 

86 to explain the observations. 

2. Observations and Theoretical Models 
2.1. Observations 

s? Bowles et al. [1963] used the radar backscattering experiment at 50 MHz at Jicamarca, 

as Peru to study the echoes of the field-aligned irregularities in the equatorial electrojet 
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89 (EEJ). Later, the experiment was repeated for different case studies by Cohen et al. [1967], 

9 0 Balsley [1969], Balsley and Farley [1971, 1973], Farley and Balsley [1973], and Fejer et al. 

91 [1975, 1976] in order to examine the characteristics of the backscattered echoes with the 

92 goal of understanding the mechanisms of the E-region plasma irregularities at the equator 

93 and high- latitudes. 

9 4 When the electron drift velocity in the equatorial electrojet exceeds the ion-acoustic 

95 speed, an echo spectrum appears simultaneously at all elevations and its Doppler shift 

96 always equals the ion-acoustic speed to a good approximation. The spectrum of the 

97 echoes is characterized by a spectral width narrower than the echoes' Doppler shift. The 

98 Doppler shift of the echoes is isotropic in the angle between the horizontal plane and the 

99 wave propagation wavevector, 6. The instability with these characteristics is called Type-I 
wo instability [Bowles et al. 1963, Balsley 1969, Fejer et al. 1975]. 

101 However, as the electron drift speed falls below the ion acoustic speed, the backscattered 

102 radar echoes exhibit different spectral characteristics and we have now what is called 

103 Type-II instability. Despite the increase of the spectral width of Type-II echoes with the 

104 wavenumber, it does not change with the elevation angle. The Doppler shift of Type-II 

105 instability varies as the cosine of the angle between the horizontal plane and the wave 
wo propagation wavevector, 9, and varies linearly with the wavenumber at constant angle 9 
10? [Bowles et al. 1963, Balsley 1969, Fejer et al. 1975, Hanuise and Chrochet 1981]. Type-II 

108 instability is always easily excited in the presence of a density-gradient [Simon 1963, Hoh 

109 1963, Maeda et al. 1963, Knox 1964, Reid 1968, Rogister and D'Angelo 1970, Balsley 
no and Farely 1971] and owing to the absence of a real electron drift threshold [Farely and 
in Balsley 1973], type-II instability is difficult to observe in the presence of type-I instability 
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112 which is excited in a strong turbulence mechanism [Sudan et al. 1983]. 

in The echoes for Type-I and Type-II instabilities can be seen both during the daytime 

ih and nighttime. The strong damping of the instabilities above 115 km, as a result of 

us the carrier's recombinations processes and the negative density scale-length, limits these 

He echoes in the daytime to the region between 93 and 113 km. However, the irregularity 

n7 in the electron density profile during the nighttime, which gives rise to irregularity in 

us the density gradient, extends the instabilities beyond 113 km and the echoes can be seen 

us up to 130 km [Kudeki 1987]. Shortly after sunset the radar backscattered data of the 

120 electrojet shows an inverse of the sign of its phase velocity. The change of the sign of 

121 the phase velocity is considered an evidence of the reversal of the equatorial electrojet, 

122 where the electrons reverse the daytime westward flow to an eastward flow during the 

123 nighttime. This electrojet reversal is accompanied by a disappearance of the instabilities 

124 echoes [Fejer et al. 1975]. 

125 Farley and Balsley [1971] shows the radar data at 50 MHz and 146-MHz frequencies that 

126 emphasize the dominance of the Type-I instability in the spectrum of the radar echoes at 
12? both frequencies when the electrons drift reaches the threshold of Type-I instability as it is 
us exceeding the ion-acoustic speed {v E > C s ). There is no, however, practical threshold for 

129 Type-II instability, as we discussed earlier, that can be detected at the HF radar echoes 

130 (50-MHz) but type-II echoes can not be seen at the VHF radar echoes (150-MHz). 

131 These measurements indicate the presence of irregularities of different scale lengths that 

132 range from a kilometer or tens of meters scales due to pure type-II instability down to a 

133 meter or shorter due to pure type-I instability. 

134 Pfaff et al. [1987] showed the data collected from sounding rocket launched from Punta, 
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135 Lobos, Peru in order to quantify the plasma instabilities by detecting the turbulence in the 

we daytime plasma density and electrostatic field in the equatorial electrojet. A simultaneous 

is? measurement taken by the Jicamarca radar showed strong type-I echoes of 3 m wavelength 

us accompanied with long-scale horizontally propagating waves. The data showed instability 

139 in the plasma density in the region of the positive vertical gradient between 90-110 km 

140 altitude with maximum variations of spectral energy measured at frequencies below 100 
ui Hz. 

«2 Later, Pfaff et al. [1997] described the measurements gathered from the in-situ sounding 

us rocket that was launched from Alantara, Brazil in 1994. Pfaff found a maximum daytime 

144 vertical electric field of about 9-10 mV/m near 105 km altitude coincident with a maximum 

us current density of magnitude 8.0fiA/m 2 . The maximum vertical electric field corresponds 

M 6 to 360-400 m/s westward electron drift, which is about the ion-acoustic speed calculated 

147 at the same altitude. These measurements agree with the estimated large electrojet Hall 

us current and Cowling conductivity estimated at about the same altitude, 

us Recently, a combined system of five radars was used at Jicamarca on July 26, 2006 

iso to monitor the equatorial electrojet, including range-time-intensity (RTI) mapping, radar 

151 imaging, radar oblique scattering, Faraday rotation, and multiple frequency scattering 

152 using AMISR prototype UHF radar [Hysell et al. 2007] . Radar imaging data shows narrow 

153 type-I echoes excited from the vertically polarized electric field of the large-scale waves 

154 which is different from the type-I echoes excited directly by the background electrojet 

155 current. Hysell et al. [2007] found stronger echoes come from the upward flux at the 
iso spectrograms from the AMISR prototype at Jicamarca and considered that as an evidence 
157 of " updown" asymmetry in the vertical particle fluxes. In addition, an asymmetry in the 
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us irregularity drifts is found in "east-west" direction with more westward drift during the 

159 daytime. 

2.2. Theoretical Interpretations 

wo Sudan et al. [1973] uses the two-step mechanism to explain the generation of the meter 

lei scale Type-II irregularities detected by 50-MHz incoherent scattering radar at electron 

162 drifts of order 100 m/s. The gradient-drift instability is excited in the daytime with 

W3 positive gradient and westward drift of the electron. When the amplitude of large-scale 

164 instability reaches a certain level, the energy starts to transfer from that primary long- 

W5 wavelength waves to secondary short-wavelength waves. This energy transfer takes place 

we due to the strong perturbations of the electrojet local parameters, where the horizontal 

167 density gradient becomes greater than the background vertical gradients and the magni- 

168 tude of the vertically perturbed drifts increases to the order of magnitude of the horizontal 
leg electron drift [Farely and Balsley 1973]. However, the excitation of vertical irregularities 
™ of wavelength of the order of meters or sub-meter is a pure type-I nonlinear mechanism 

171 that takes place only when the electron drifts exceed the ion-acoustic speed [Sudan et al. 

172 1973]. 

173 On the other hand, the dependence of the phase velocity of type-I instability on the 

174 cosine of the angle between the horizontal plane and the wave propagation wavevector, 

175 9, in the linear theory is not supported by the observations [Sudan 1983a]. The phase 

176 velocity which equals the Doppler shift is independent of the angle 9 and has a constant 

177 value that equals the ion-acoustic speed [Balsley 1969, Farely and Balsley 1973, Hanuise 

178 and Crochet 1981]. Kamanetskaya [1971] and Rogister [1971], Kaw [1972], Skadron and 

179 Weinstock [1969], Weinstock and Sleeper [1972], Rogister and Jamin [1975], and Jamin 
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wo and Kennel [1976] suggested various mechanism to explain these observations such as the 

181 quasilinear theory that limits ve to C s , the refraction of waves in the electrojet velocity 

182 profile, the nonlinear orbital diffusion of ions or electrons, mode coupling outside the in- 
ns stability cone, and the propagation with finite k\\ along the magnetic field, respectively, 
i8 4 but they all rely on incorrect hypotheses by either treating type-I and type-II instabilities 
las independently [Sudan 1983a] or considering homogeneous electrojet [Farely and Balsley 
ise 1973]. 

is? Sudan [1983a, 1983b] has also treated this problem differently by considering the co- 

188 existence of the two types of instabilities and the effect of the developed type-II strong 

189 turbulence on the evolution and development of type-I instability. Sudan [1983b] sug- 
wo gested that the nonlinear process modifies the electron-neutral collision frequency u en by 

191 wave-amplitude dependent value that causes the growth rate to vanish in the saturation 

192 region. 

3. Plasma Dynamics and Model Equations 

193 The E-region is bounded by two layers of reduced conductivities; below 90 km the 

194 collision frequencies of the charge carriers with background neutrals are very high, reducing 

195 the conductivity to very small values, whereas above 130 km both the electrons and 
we the ions are magnetized and are influenced by the Lorentz force which also reduces the 
197 conductivity to very small value. According to the sounding rocket data, the maximum 
us conductivity is centered at 105 km. Hence, any vertical current will be inhibited by a 

199 vertically induced electric field as a result of the accumulation of the electrons at the 

200 upper layer and the accumulation of the heavy ions at the lower layer. We adopt a slab 

201 geometry model with the x-axis pointing northward along the geomagnetic field which is 
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202 perpendicular to the plane in which the dynamical equations are written, y-axis pointing 

203 westward, and the z-axis pointing vertically upward, as shown in figure(l). 

204 The possible nonlinear coupling mechanisms between the linear unstable modes is 

205 shown in figure (l-a,b). The forward-cascade coupling mechanism in figure (1-a) shows 

206 the interaction between nearly collinear modes to produce a mode of shorter wavelength. 

207 However, the interaction between two unstable modes of short wavelength to produce a 

208 stable mode of longer wavelength in an inverse-cascade coupling mechanism is shown in 

209 figure (1-b). Kudeki and Farley [1989] and Lu et al. discussed the sensitivity of the aspect 

210 angel of the equatorial electrojet irregularities and the coupling mechanism between the 

211 linearly unstable modes. 

3.1. Dynamic System Equations 

In our model, we have three dynamical equations for three fields; the plasma density 
n, the electric potential <p, and the ion velocity potential x- The ion density profile of all 
species shows that is the dominant ion species in the E-region, and we consider only 
one ion species. We assume isothermal state, T e = Tj, and quasineutrality, n e = n«, since 
the electron Debye length, Ad ~ 1cm, is much smaller than the spatial variation of the 
electron density. The small ratio between the ion gyrofrequency, u C i, and the ion-neutral 
collision frequency, is in , makes the mean-free path of the ions small so that ions cannot 
complete one gyrocycle before experiencing another collision. Therefore, the ions are 
considered unmagnetized and we ignore the Lorentz force, v\ x B; in the ion's equation 
of motion. However, the electrons are considered magnetized since the ratio between 
the electron gyrofrequency, u ce , and the electron- neutral collision frequency, v in , is much 
larger than unity [Kelley 2009] . Finally, the acceleration due to gravity in both equations 
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of motion is ignored. 

The negligible variations of the equatorial geomagnetic field, 8B 0 = 20nT, with respect 
to the background geomagnetic field, B Q = 30/iT, makes us consider an electrostatic 
plasma with the electric field is defined as E = — V</?. Also, the ion fluid velocity can be 
written as Vi = — 

The ion's continuity equation and the equations of motion of electrons and ions can be 
written in their general forms as: 

d t n + V • (nvi) = 0 (1) 

rriinAtVi = eriE — TjVn — V ■ 7r — minu in Vi (2) 

m e nA t v e = — en(E + v e x B) — T e Vn — m e nu en v e (3) 

where, A t is the convective derivative; A t = (<9 t + v • V), and fr is the ion viscosity stress 
tensor and its elements are given by: 



nTi 

n a p = 



(dpv a + d a v/3 - jjv • v5 a ^j (4) 



In this model we consider the ion viscosity and the electron inertia. The physics of these 
two terms and their effects on the stabilization of the growth rate and the saturation of 
the instabilities will be explained later. 

Equations (1) and (2) can be rewritten in terms of the normalized quantities (n = n/n 0 , 
<p = e<p/T e , and x = x/ 1 ) as: 



d t h = Vw • V% + nV x (5) 
d t V 2 x = V 2 (y\h + v\q> - v m x + (Vx) 2 ) + ^V 4 x (6) 
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212 The magnetized electron drift velocity, in a plane perpendicular to the geomagnetic 

213 field, can be given after some algebraic manipulation of equation(3) as 

214 

if — ^finn , 

and for the first-order approximation we can rewrite the electron velocity as 

v e ± =v E + v de + v Pe + v v (7) 

216 where, ve is the E x B drift velocity, Vd e is the diamagnetic drift velocity, vp e is the 

217 polarization drift velocity, and v u is the drift velocity due to the frictional force between 
2ia the electrons and the neutral background. 

219 Substituting the ion drift velocity, Vi = — Vx, and the electrons drift velocity from 

220 equations (7) in the plasma quasineutrality condition, V-(Ji + J e ) = 0, and by making 

221 some algebraic manipulation along with using the normalized quantities mentioned above, 

222 we can close our system with the following normalized dynamic equation: 

d t V 2 <p = v en Vh ■ Vn - v en Vh ■ + p^ 2 Vn • V% (8) 

+z/ en V 2 n - u en V 2 (p - p- 2 V 2 x 

223 where, [/, g] is the Poisson (Jacobian) bracket which is defined as [/, g] = d x fd y g — d x gd y f. 

3.2. Instabilities Evolution and Saturation 

224 The ion viscosity term in equation(6), f ^-Vx, plays an important role in stabilizing 

225 the modes at the large k-values and it simulates the effect of Landau damping, which can 

226 not represented in the fluid models. The effect of this term on the growth-rate is discussed 

227 in more detail in section(4). 
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228 In equation(8), the Jacobian — Q ce 5ip, 5n has a coefficient that is two- and four-orders 

229 of magnitude larger than the coefficients of the other terms, and it shows the coupling 

230 between the horizontal and the vertical waves in the electrojet. In addition, this Jaco- 

231 bian bracket is responsible for the generation of the vertical structures with the down- 

232 ward/upward drifts of the low/high density plasma which is driven by the gradient-drift 

233 instability. 

234 On the other hand, the Jacobian bracket —p 2 e Vl ce 5(p, V 2 <5y? , which considers the ExB 

235 drift of the plasma vortices, adds more damping to the growing modes at high k-values 

236 to drive the system to saturation when the electron density amplitude grows linearly to 

237 (12%-25%) of its initial value at the beginning of the simulation, depending on the solar 

238 conditions. 

239 Hence, the dynamics of the electron density, electric potential, and the ion velocity po- 

240 tential of the ionosphere plasma at a plane perpendicular to the geomagnetic field in the 

241 equatorial E-region can be described by a set of hyperbolic partial differential equations 

242 (5), (6), and (8). 

4. Linear Regime 

243 First we study the basic properties of the instabilities in the linear regime. 

4.1. Dynamic System Linearization 

To study the linear characteristics of the dynamical system, n, ip, x, we need to linearize 
the equations (5,6,8). The linearized quantities for density n, electric potential <p, and 
ions drifts potential Xi with the assumption of a stationary background plasma Vxo = 0, 
are: 
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n = n Q + Sn, (p = (p Q + 5(p, \ = $X 

where the general spectral representation of the fluctuations of the perturbed part of those 
three fields in space and time is given by: 

5<p(i = 1, 2, 3) = J2k x ,k y ^ exp [i(k • x - ut] 
with a condition for real values giving <p n (—k, —u*) = (p* n (k,uj). where (501,502,503) are 
equivalent to (5n,5ip,Sx), respectively. 

After dropping the zero-order and higher order terms the set of linearized and normalized 
equations is: 

d t 5h = V 2 5x (9) 

244 

d t V 2 5<p = (10) 

/ 2 v E v E v en 2z/ en \ 
v en ^ -d y + — d z + -j—d z Sn 

~ UnV 2 + V E d y V 2 + ^d y + ^d z ) Sip 

d t 5x = v 2 M + v\b(p - L n - ^V 2 ) 5x (11) 

245 where p e = v e /u ce is the electron gyroradius, L n = d z £nn 0 is the electron density scale- 

246 length that has an approximate value of 6 km (1 km) during the daytime (nighttime) with 

247 n 0 as the background average charge- carrier density, and v E = —B~ 1 d z (p 0 is the Ex B 

248 drift speed 400.0 m/s) with (p„ as the background average electric potential. 

249 The analytical solution of the set of three linearized equations of the dynamical system 

250 for the three normalized quantities gives rise to a cubic dispersion equation as shown in 
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25i equation(ll): 



4.2. Numerical Solution of Linear System 

252 The local and the nonlocal approximations have been used to model plasma instabilities 

253 in the E-region depending on the relation between the density scale-length, L n and the 

254 fluctuation wavenumber, k [Ronchi 1989]. We use the local approximation, valid in the 

255 case of \kL n \ 3> 1 (or L n 3> A, where k = 2ir/\), which is the condition that is satisfied 

256 in this study for the length of the type-II irregularities smaller than a kilometer. 

257 We employ the the International Reference Ionosphere (IRI2007), Naval Research Lab- 

258 oratory Mass Spectrometer and Incoherent Scatter Radar Exosphere (NRLMSISOO), and 

259 International Geomagnetic Reference Field (IGRF11) empirical models under different 

260 conditions to get the altitude-dependent values of the collision frequencies, gyrofrequen- 

261 cies, and gyroradii for ions and electrons, as shown in table(l) for the solar-maximum and 

262 solar-minimum conditions. 

4.3. Altitude Dependent Growth Rate 

263 To study the altitude-dependence of the instabilities in our unified model we solve 

264 numerically the eigenvalue problem of the set of linearized equations (9,10,11) between 

265 90-120 km to find the growth rates and phase velocities at each altitude. We use 0.03 

266 rrr 1 resolution in k-space and 0.25 km resolution along the altitude, z-axis. 

267 The growth rate profile as a function of altitude is shown in figure(2) and divides the 




(12) 
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268 E-region into three regions. The lower region below 103 km shows a dominance of the 

269 pure gradient-drift (type-II) instability over the small wavenumbers with the absence of 

270 Farely-Buneman instability at the intermediate and large wavenumbers due to the large 

271 collision rates of plasma species with the ionosphere neutral background and the small 

272 relative speed between the electrons and ions. The region extends between 110-115 km 

273 does not have instabilities at small wavenumbers which reflects the absence of gradient- 

274 drift instability as a result of the negative scale-length of the plasma density. The growing 

275 modes in this region are found in the intermediate and large wavenumbers where the pure 

276 Farely-Buneman (type-I) instability dominates. A coupling between type-I and type-II 

277 instabilities can be seen at the core of the electrojet region between 103-110 km, where 
27a the maximum positive growth rate extends over all wavelength scales from hundreds or 

279 tens of meters to sub-meter. 

280 Thus, the altitude-dependent growth rate profile is in agreement with the rocket obser- 

281 vations that divides the instabilities in the equatorial electrojet region into three regions 

282 as we described with the peak of the electrojet turbulence between 103-108 km [Pffaf 

283 1987b]. 

4.4. Standard Model versus Unified Model 

284 In figure(3) we demonstrate the effect of the ion viscosity term, V 4 x, and the electron 

285 inertia term, d y V 2 (p, which both are stabilizing at high wavenumbers for a reference case 

286 at 105 km altitude. When the electron inertia is ignored but the ion viscosity is retained 

287 (dashed-red) , we find the instabilities are suppressed at high wavenumber (k y ~ 25m" 1 ). 

288 However, the stabilization of the growing modes is found at k y ps 81m" 1 when we only 

289 consider the electron inertia term and drop the term coming from the ion viscosity. The 
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290 electron inertia is not relevant at these length-scales in the linear regime, however, it is 

291 important in saturating the growing modes of large amplitudes in the nonlinear regime. 

292 A monotonic increase of the growth rate with the wavenumber (dashed-green) can be 

293 seen when we drop both the ion viscosity and the electron inertia terms. This case is 

294 similar to the standard two-stream fluid model [Fejer et al. 1975] that shows no cut- 

295 off for the type-I instability. However, the solid black line shows stabilization at large 

296 wavenumbers around k y = 11.5(m _1 ) when both the ion viscosity and electron inertia are 

297 included in the model. 

298 A maximum growing mode at k y =7.8(m -1 ) is equivalent to an irregularity of scale size 

299 less than a meter in the zonal direction. In the vertical direction we have the maximum 

300 growing mode at k z =0, and there are modes of positive growth rate at k z =2.1(m~ 1 ), 

301 which is equivalent to 3 meters wavelength. Farely and Balsley [1973] found that the 

302 excitation of 3 meters wavelength waves requires a very large drift speed which is not 

303 possible in the equatorial electrojet region and they referred to Sudan's et al. [1973] 

304 theory, which attributed the generation of vertical structures of wavelength 3-meters or 

305 less in the equatorial electrojet to the cascade of energy to small scales. 

306 The linear growth-rate for both types of instabilities that is derived in the standard fluid 

307 model by Fejer [1975] increase monotonically in contrast with the growth rate calculated 

308 numerically from our unified fluid model that includes both the ion viscosity and electron 

309 inertia terms. Our unified model produces a cut-off at high wavenumbers. The growth 

310 rate profile of the unified fluid model is similar to that of the kinetic model of type-I 

311 instability studied by Schmidt and Gary [1973] and Oppenheim et. al. [1996], which was 

312 compared to the linear growth rate of the standard fluid model by Fejer [1975]. 
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so On the other hand, there is a decrease in the phase velocity found in our unified model 

314 compared to that which was estimated by the standard linear fluid model [Fejer 1975] 

315 as shown in figure (4). This decrease in the phase velocity at high wavenumbers is in 
3W agreement with the radar observation that shows a phase velocity limited to the acoustic 
si? speed of the ions inside the equatorial electrojet. 

4.5. Type-I and Type-II Instabilities 

sis The linear growth rate profile shows the presence of a transition from Type-II insta- 

3w bility at small wavenumber (< 0.1m -1 ) to the Type-I instability at large wavenumbers. 

320 For electron drift velocity less than the ion-acoustic speed we find that the growth rate 

321 maximum is located at the transition wavenumber 0.1m -1 ) with a peak value that 

322 depends on the electron density scale length (L n ). Also, the level of the growth rate 

323 peak at the large wavenumbers, where type-I instability dominates, depends on both the 

324 electron drift speed (ve) and the electron density scale length (L n ). Therefore, to include 

325 both instabilities, type-I and type-II, in our simulation we need both a small minimum 

326 wavenumber and a large box in k-space. We discussed the system initialization in the 

327 next section. 

5. Nonlinear Regime and Results 

328 Nonlinear evolution and saturation of Type-I and Type-II instabilities have been studied 

329 for several decades. The one-dimensional simulations that have been done by Sato and 

330 Tsuda [1967] and Rogister [1972] to study the gradient-drift (Type-II) instability were 

331 criticized by Sudan [1973] because they ignored the coupling between the waves prop- 

332 agate horizontally and vertically in the equatorial electrojet region. The saturation of 



DRAFT 



January 14, 2015, 12:23pm 



DRAFT 



X - 20 EHAB HASSAN ET AL.: EQUATORIAL ELECTROJET TURBULENCE: BASELINE2D 

333 Farley-Buneman (Type-I) instability has also been studied in different types of simula- 

334 tions. Newman and Ott [1981] ran a fluid nonlinear simulation that shows the dominant 

335 wave to propagate in a direction that is different from the direction of the electron drift. A 

336 particle-in-cell (PIC) code was used by Machida and Goertz [1988] to model the dominant 

337 wave in a plane parallel to the geomagnetic field. This geometry ignores the nonlinear 

338 term, b ■ Vy? x Vw where b is in the direction of the magnetic field, that is important for 

339 the saturation of Type-I instability. 

3M Oppenheim et al [1996a] presented two-dimensional hybrid simulations that treat the 

341 electrons in the fluid regime and the ions kinetically to model the effects of Landau damp- 

342 ing. The agreement between the radar/rocket spectra and the simulation spectra was 

343 presented by Oppenheim and Otani [1996b]. The hybrid simulation is able to reproduce 

344 the phase velocity during the saturated state similar to the drift velocity of the irregulari- 

345 ties and its independence from the elevation angle. Also, the hybrid simulation shows the 

346 presence of a secondary Type-II instability in the absence of the electron density gradient 

347 as a result of the mode coupling. Otani and Oppenheim presented the mechanism of the 

348 saturation of Farley-Buneman (Type-I) instability in a series of papers [Otani and Oppen- 

349 heim 1998, Otani and Oppenheim 2006]. Otani and Oppenheim attribute the saturation 

350 mechanism to the interaction between three modes in the system, and used a fluid model 

351 to reproduce some observed features from the echoes of the irregularities. 

352 In 2008, Oppenheim et al. presented the results of a large-scale simulation from a 

353 fully kinetic simulator called Electrostatic Parallel PIC or EPPIC code. This is a high- 

354 resolution two-dimension simulation with (4096 x 4096) grid points that covers (160 m 

355 x 160 m) space which is comparable to our simulation box. The simulation showed the 
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356 coupling between the linearly growing modes and the damped modes and included the 

357 electron and ion thermal effects. This simulation reproduced most of the measurements 
sss and features of Farley-Buneman turbulence with a phase velocity larger than the ion- 

359 acoustic speed. Later, Oppenheim and Dimant [2013] presented the first fully kinetic 3-D 

360 simulation results for the Farley-Buneman instability with the ionospheric condition at 

361 the aurora region. The simulator uses a sophisticated method in calculating the collision 

362 frequency of the ions and electrons, in addition to the anomalous heating of the electrons. 

363 The simulation was able to resolve the 1-5 meters waves during the transition region and 

364 the formation of a large structure at the saturation region. Also, it shows the role of the 

365 anomalous electron heating in reducing the phase velocity of the short waves to a value 

366 very close to the local ion-acoustic speed compared to the previous 2-D simulation. The 

367 electrons and ions thermal effects are out-of-scope of the current work in this paper and 
sea they will be considered in future work. 

5.1. Nonlinear System Initialization 

369 In our nonlinear simulation, we use the pseudospectral method to find the spatial deriva- 

370 tives of the fields in k-space and then integrate the three differential equations in the time 

371 domain using a sixth-order Runge-Kutta (RK6) method with periodic boundary condi- 

372 tions at all the boundaries of all fields. We solve our nonlinear system in a simulation box 

373 of size 100 m x 100 m in the r-space with grid resolution 1024x1024. This r-space grid 

374 is equivalent to 22 m _1 x 22 m _1 in k-space (after removing the dealiasing region that 

375 is required by the nonlinear terms for spectral methods). We use 1024 x 1024 resolution 

376 to be able to include both types of instabilities (Type-II at k y ^0.1 m~ l and Type-I at 

377 k y ^7.8 m" 1 ) in the simulation box of our dynamical system as we indicated earlier in 
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378 section (4) and show in figure (3). 

379 The local values (where, \kL n \ 3> 1) around 105 km altitude for the physical quantities 

380 on March 2008; which is coincident with the deep minimum and a solar quiet (Sq) time of 

381 the last and unusual solar minimum conditions [Benevolenskaya and Kostuchenko 2013], 

382 with setting ve = 400m/s and L n = 6km are employed in the coefficients of the dynami- 

383 cal system equations. Woods et al. [2010] found the total solar irradiance in 2008 is much 

38 4 smaller than its value in the previous solar minimum in 1996, and the value of the solar 

385 radio flux index at this day is F10. 7=69.9. 

sse We initialize the plasma density, electric potential, and the ion velocity potential in 

387 k-space using the values of the eigenvectors calculated numerically in the eigenvalue prob- 

388 lem at all the k y and k z values corresponding to positive growth rates as discussed in 

389 section(4), and we run the simulation for more than 1 second with 10 /xs time step. 

5.2. Spectral Analysis of Electron Density 

390 In the nonlinear run we evolve the perturbed normalized quantities (5n, Sip, Sx) in k- 

391 space over time to check the transition from the region of linear dominance to the region 

392 of nonlinear dominance and end in the saturated state of the simulation. Snapshots of 

393 the energy spectrum of the plasma density are shown in figure(5) at different stages of the 

394 simulation. For context, the reader is also referred to figure (3) (which shows the range of 

395 linear instability), and figures(6-8) (which show the time evolution of the density pertur- 

396 bations). 

397 The first panel in figure (5) shows the energy spectrum of the plasma density in the 

398 linear growth phase with more energy content at small k-values in the horizontal and 

399 vertical directions. In the second panel, the energy is concentrated at low k z and shifts to 
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00 higher k y to damp linearly after generating waves of sub-meter wavelengths. Panels three 

01 through five demonstrate the transition from initial nonlinear stabilization to fully satu- 

02 rated turbulence. The third panel (t = 0.53 s) represents a point of transient saturation, 

03 as can be seen in figure(10). At this point, the energy is largely concentrated in the small 

04 type-I scale range. In panel four, the dominant scales are reversed during the transition 
os to the final saturated state. Here the energy is concentrated at large type-II scales, with 
oe comparatively less energy at small scales. In the fifth panel, the final nonlinear balance is 
07 achieved with two distinct scales having energy apportioned quite evenly between them, 
os These two dominant scales, however, are embedded in broadband turbulence, as can be 

09 seen from the significant energy observed over the entire instability range (e.g., in the 

10 intermediate k y = 2.0 — 7.0 range). 

11 The onset of the nonlinearity dominance in the transition region is coincident with the 
is generation of long structures of the plasma irregularities in the vertical direction with 
o scale length of 40 meters, see figure (7). In addition, we attribute the generation of the 
a 3 meter scale waves in the vertical direction to the cascade energy transfer from the long 
is structures at small wavenumbers to the short irregularities at the large wavenumbers 
i6 which is consistent with the theory of Sudan et al. [1973] for the generation of small-scale 
I? irregularities in the equatorial electrojet. 

is This complex nonlinear interplay between the two scales and the back and forth exchange 

19 of energy between them will be studied in detail in future work using nonlinear energy 

20 transfer functions. 

5.3. Instabilities Evolution, Coupling, and Saturation 
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421 Many of the features described in the previous section can also be directly observed in 

422 real-space contour plots. The upward (downward) drift of the low-density (high-density) 

423 plasma is shown in figure(6) during the linearly growing phase of the simulation. The ver- 

424 tical drift velocity is found from successive frames to be around 50-70 m/s. The plasma 

425 is broken horizontally into small-scale waves of 1-3 meters in wavelength. 

426 Figures(7 and 8) show snapshots of the plasma density fluctuations at t — 0.6 s (during 

427 the initial saturation phase) and t = 1.0 s (during the final saturated phase), respectively. 

428 The multi-scale nature of the turbulence is evident in the East- West direction in both fig- 

429 ures, with prominent features visible at ~ 1 m that are embedded in structures with ~ 10 

430 m scale. In the vertical direction, the fluctuations exhibit a stark change during this tran- 

431 sition to fully developed turbulence. The vertical structures are very large scale ~ 20 — 40 

432 m at t — 0.6 s, but are reduced to ~ 1 — 3 m at t — 1.0 s. At this point, the simulation 

433 box is filled with small-scales structures embedded both vertically and horizontally in the 

434 long-scale irregularities of Type-II instabilities. The small-scale vertical structures can be 

435 considered as the formation of secondary type-I instabilities in the vertical direction and 

436 is consistent with the energy cascade theory [Sudan et al. 1973] to small vertical scales. 
43? A closer look at the turbulence in the plasma density under the solar-maximum con- 

438 dition (March 1987) is shown in figure(9). The left-panel shows the long-structures at 

439 t = 0.75 s which is right-shifted in time with respect to the transition region under the 

440 solar-minimum condition as we discussed in the next section. These long-structures fill 

441 the box with a few small-scale structures. At t = 1.2 s, the perturbation in density is 

442 saturated and we see in the right-panel of figure(9) many small-structures of 1-3 meter 

443 scale-length. The structure of the plasma density during the growing-, transition-, and 
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4M saturation-phases of the simulation is identical under different solar conditions and only 

445 the saturation level of the perturbation is different. 

446 Watching successive frames of the density evolution during the saturation states shows 

447 that the downward and upward drifts are not in a straight vertical direction, but they 

448 drift with an inclination to the horizontal direction which is parallel to the direction of 

449 ExB plasma drifts. Oppenheim and Otani [1996] called this motion toward the corners of 

450 the simulation box "a bent of the primary waves," and they attributed it to the decrease 

451 in the linear growth-rate for the primary waves of k z ^ 0 to a value that is below its value 

452 at k z = 0 in the linear regime of the simulation. This feature can be seen from the vertical 

453 profile of the growth rate. 

5.4. Turbulence in Electron Density 

454 A comparison between the temporal variation of the plasma density during the condi- 

455 tion of solar maximum on March 1987 (blue-line) and the condition of solar-minimum on 

456 March 2008 (red-line) is shown in figure(lO). The top-panel shows the maximum pertur- 

457 bation in the plasma density, however the bottom-panel shows the standard-deviation of 

458 the perturbed density as a function of time. 

459 The standard-deviation shows distinctly the three domains of the simulation: linear, 

460 transition and saturation. The evolution time required in the simulations in the growing 

461 phase under the solar-maximum conditions is longer than the corresponding time under 

462 the solar-minimum conditions. Also, the level of the standard-deviation of the saturated 

463 perturbed plasma density is higher under the solar-maximum conditions (a$ n = 3.0%) 

464 than the solar- minimum conditions {a^ n = 1.7%) with a At = 0.1 s difference between 

465 their start of the saturation region in the simulation. 
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466 Under the solar-minimum condition, the perturbation in the electron density increases 

467 to 12% of the background density during the transition region before it goes back to 7% 

468 in the saturation region as a result of the strong coupling between the primary and sec- 

469 ondary waves in the nonlinear saturation regime. 

470 This is below the expected 10-12% perturbations found in both the Rocket and the 

471 Radar Observations [Balsley and Farely 1971, Kudeki et al. 1987, Pfaff et al. 1997, Hysell 

472 et al. 2007] . Our interpretation of the low level of plasma density perturbations is the re- 

473 suit of the low solar radio flux accompanied by the unusual solar minimum conditions and 

474 the solar quiet condition employed in this simulation [Benevolenskaya and Kostuchenko 

47 5 2013, Woods et al. 2010]. 

476 When we examined the dynamic system under solar-maximum conditions we found 

477 that the temporal variation of the perturbed density during the saturation region shows 

478 a maximum much higher than the 7% found during the solar-minimum conditions. The 

479 perturbation in the electron density increases to 25% of the background density during the 

480 transition region before it goes back to 15% in the saturation phase. In a following paper, 

481 we are going to examine the effect of the background electric field and density-gradient 

482 scale-length on the maximum value of the perturbed density. 

5.5. Perturbed Electric Fields 

483 The maximum perturbed electric field in the zonal and vertical directions are shown 

484 in the upper and lower panels of figure(ll), respectively. The zonal component of the 

485 perturbed electric field grows to 140 mV/m before it saturates at 80 mV/m under the 

486 solar-maximum conditions which is double its saturation value under the solar-minimum 

487 conditions. Similarly, the vertical electric field saturates at a larger value (35 mV/m) 
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488 compared to the solar minimum saturation value (15 mV/m). The difference in the per- 

«9 turbed field between the solar maximum and solar minimum is directly related to the 

490 level of the perturbed density under these different conditions which is discussed in the 

491 previous subsection. 

492 The root-mean-squares (rms), which are measurable quantities, for the perturbed elec- 

493 trie field in the zonal (solid-lines in upper-panel) and vertical (solid-lines in lower-panel) 

494 directions are shown in figure(12). While the magnitude of the rms of the zonal perturbed 

495 electric field in the saturation regime reaches 15 mV/m during the solar maximum con- 

496 ditions, it does not exceed 8 mV/m under the solar minimum conditions. On the other 

497 hand, the rms value (5E Z = 3.0 mV/m) of the perturbed vertical electric field under the 

498 solar-maximum conditions is almost double its value (5E Z = 1.7 mV/m) under the solar- 

499 minimum conditions. 

boo The simulation results of the magnitude of the zonal and vertical perturbed fields under 

sol the ionosphere conditions of March 1987 shows a good agreement to the measured rocket 

502 data in the CONDOR campaign of the perturbed electric field in the zonal direction. 

5.6. Plasma Fluxes Asymmetries 

503 The "east-west" zonal and "up-down" vertical components of the perturbed electric 

504 field are shown in the upper and lower panels of figure(12), respectively, during the solar- 
sos maximum (red) and solar-minimum (blue) conditions. While the solid lines represent 
506 the rms of total perturbed electric fields, the dotted and dashed lines represent the rms 
so? of the westward (upward) and eastward (downward) perturbed components in the zonal 

508 (vertical) directions, respectively. 

509 These zonal and vertical components of the electric field give rise to Pedersen drifts 
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510 along its direction and Hall drifts in a direction perpendicular to the plane containing 

511 both the electric field and the magnetic field as a result of the ExB drift. The large ratio 

512 between the Hall and Pedersen conductivities makes the effect of the Hall current more 
so important in the equatorial electrojet region. In turn, the difference in the Hall currents 
5« in the zonal and vertical directions gives rise to an asymmetry in the drift velocities and 
sis density fluxes in both directions. 

sie The large difference between the upward and the downward components of the per- 

517 turbed electric field is common in all phases of simulation and causes an "east- west" 

sis asymmetry in the drift velocity. Patra et al. [2005] reported the first observations veri- 

519 fying the "east-west" velocity asymmetry found in Type-II echoes using a 18 MHz radar 

520 located at Thumba Equatorial Rocket launching Station (TERLS), Trivandrum, India. 

521 Patra et al. [2005] attributed the "east- west" velocity asymmetry to the large-scale pri- 

522 mary waves that give rise to a large upward plasma drift. However, we find the "east- west" 

523 asymmetry exists even after the breaking-up of the large-scale vertical waves into smaller- 

524 scale structures. Thus, we consider the "up-down" asymmetry of vertically perturbed 

525 electric field, which is responsible of the zonal Hall drifts, to be the cause of the "east- 

526 west" velocity asymmetry. 

527 On the other hand, the asymmetry in the vertical fluxes produces echoes of different 

528 strengths in the spectrograms from the AMISR prototype at Jicaqmarca as reported by 

5 29 Hysell et al. [2007] and studied by Fejer et al. [1976], Farley et al. [1978], and Kudeki et 

530 al. [1985]. The Hall drifts in the vertical direction do not show large asymmetry because 

531 of the very small difference between the zonal components of the perturbed field as shown 

532 in the upper panel of figure(12). The large difference between the zonal components can 
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533 only be found in the presence of large-scale structures during the transition phase of the 

534 simulation, under the condition of solar maximum, and then decreases as it goes deep in 

535 the saturation phase with the "break-up" of the large-scale irregularities into small-scale 
sse structures. However, the large difference in the Pedersen drifts may play a role in causing 
537 the asymmetry in the density flux along the vertical direction. 

5.7. Phase Velocity of Plasma Irregularities 

sss The plasma density in the zonal direction as a function of time is shown in figure (13). 

539 The inverse Fourier transform of the density spectrum is taken for all k y values of k z =0 

540 (m _1 ) for all time steps in the saturation state of the simulation. We can see small-scale 

541 waves of Type-I superimposed on the large-scale waves of Type-II. The phase velocity 

542 of the small-scale waves has a value of 330 (m/s) and drifts westward. This value of 

543 phase velocity is very close to the ion-acoustic speed (318-323 m/s) calculated from the 

544 ionosphere background data. So, the phase velocity of the irregularity waves is limited to 

545 the ion-acoustic speed at the core of the electrojet. This result shows a good agreement 

546 with the radar and rocket observations that limit the speed of the plasma irregularities in 

547 the ionosphere E-region to the acoustic speed of ions and show them propagating in the 

548 westward direction [Farely and Balsley 1973, Kudeki et al. 1987]. 

6. Summary and Conclusion 

549 The temporal and spatial evolution of the plasma density, electric potential, and ion 

550 velocity potential are studied linearly and nonlinearly in a two-dimensional fluid model. 

551 Our unified fluid model considers both ion viscosity and electron inertia. The dominant 

552 nonlinear terms ,[5(p, Sn] and [dip, V 2 50\, are responsible for the coupling between Type-I 
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553 and Type-II fluctuations and drive the dynamic system to the statistically steady state. 

554 In our unified model, the electron inertia, which has been ignored in the earlier mod- 

555 els, adds the polarization drift to the dynamic system and is important for the nonlinear 
sse coupling between different modes and the saturation of growing modes of small-scale 
ss? structures. The ion viscosity models Landau damping, which, in contrast with kinetic 
ssa models, is not included in fluid treatments. 

559 The linear results also show that the ion viscosity and electron inertia are important for 

seo stabilizing the large wavenumbers around 3 meter scale sizes in the vertical direction and 

561 less than a meter in the horizontal direction. The unified fluid model linear growth rate 

562 is comparable to the growth rate calculated with a kinetic model by Gary and Schmidt 

563 [1973] and a hybrid model by Oppenheim et al. [1996] for the Farley-Buneman (Type-I) 

564 instability. The linear phase velocity predicted from our unified model is smaller than 
ses that obtained from the standard two-stream fluid model in such a way that is generally 
see consistent with the observed values of the phase velocity and close to the ion-acoustic 
se? speed. 

sea The altitude-dependence of the growth rate profile shows the role of Type-I and Type-II 

569 instabilities in exciting the irregularities in the E-region between 90-120 km as suggested 

570 from the rocket data studied by Pffaf [1987]. Type-II instabilities dominate between 90- 

571 103 km where the plasma is very collisional and electron drifts are too small to excite 

572 Type-I instability. However, Type-I instability dominates between 110-115 km in a region 

573 of negative gradients that suppresses Type-II instability. The coupling between Type-I 

574 and Type-II instabilities takes place at the core of the electrojet region between 103-110 

575 km where the maximum growth rate is extending over a wide range of wavenumbers. 
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576 The time-series of the spectrum of the density perturbation shows the back and forth ex- 

577 change of energy between plasma irregularities of different scales. This exchange of energy 

578 can be considered to be a verification of the Sudan et al. [1973, 1983] theory of cascade 

579 energy transfer to structures of 3 meters scale in the vertical direction. The perturbed 
Sao density spectrum exhibits the coupling between Type-I and Type-II instabilities and shows 
sal the generation of long-scale structures in the vertical direction. It also demonstrates how 

582 the energy is transferred back from the large-scale perturbed density structures excited 

583 by type-II modes to the small-scale irregularities exhibiting type-I features. The plasma 

584 density spectrum shows an even distribution of energy between the large-scale structures 
sas (small wavenumbers) that produced type-I modes and the small-scale structures (large 
sse wavenumbers) that produced by type-II modes. The energy content of the irregularities 

587 of scale-length less than a meter in the vertical direction is small, which explains the diffi- 

588 culties to detect them by coherent scattered radars [Farley and Balsley 1973]. In addition, 
sag the embedding of the small-scale structures in the large-scale structures demonstrates the 

590 difficulty to observe the long-scale structures of type-II while type-I instabilities are ex- 

591 cited. 

592 Our nonlinear model is able to reproduce a number of features of rocket measurements 

593 and radar observations under different solar conditions such as: (1) a saturation of the 

594 density perturbation between 7-15%, (2) root-mean-square values of the horizontal com- 

595 ponent (5E y = 8.7-15 (mV/m)) and vertical component (SE y = 1.7-3 (mV/m)) of the 

596 perturbed electric field, (3) a linear and nonlinear reduction of the phase velocity of the 

597 horizontal westward irregularities to values less than or equal the ion-acoustic speed, (4) 

598 an asymmetry in the plasma fluxes in the vertical direction as a result of an asymmetry 
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599 of the perturbed zonal electric field components, (5) an asymmetry in the plasma drifts in 

600 the horizontal direction as a result of an asymmetry of the perturbed vertical electric field 

601 components, (6) a break-up of the large-scale vertical plasma structures into small-scale 

602 structures of 3 meters and less in the saturation state, (7) an inclination of the small 

603 structures drift toward the corners of the simulation box as a result of the vertical up 

60 4 (down) drifts of the depletion (enhanced) density regions due to Type-II instability. 

605 The low-level of fluctuations of the electron density and the small eastward zonal per- 

606 turbed electric field are attributed to the solar quiet (Sq) and unusual minimum conditions 

607 with low total solar irradiance on March 2008. 

60s In conclusion, this unified fluid model is able to reproduce comparable results in the 

609 linear regime to the kinetic theory for both the growth rate stabilization and the reduction 

610 in phase velocity. It also exhibits quantitatively the coupling between Type-I and Type-II 
en instabilities excited in the equatorial electrojet region and the cascade transfer of energy 

612 to small structures as proposed by Sudan [1973]. The excitation of secondary Type-I 

613 instabilities from Type-II is verified and seen from the spectrum of the perturbed density 

614 at 105 km altitude. Finally, the nonlinear results match most of the radar observations 

615 and rocket measurements. 

616 Employing an advanced parallel computing technique such as Message Passing Inter- 
in face (MPI) with this stable fluid model will ensure the retrieve of all the E-region features 
6i8 with a large-scale simulation that covers equatorial electrojet layer. 

6w Acknowledgments. This work is supported by the National Science Foundation grant 

620 AGS-0964692 to the University of Texas at Austin. The first author has been partially 



DRAFT 



January 14, 2015, 12:23pm 



DRAFT 



621 
622 



EHAB HASSAN ET AL.: EQUATORIAL ELECTROJET TURBULENCE: BASELINE2D X - 33 

supported by Ain Shams University, Cairo, Egypt and the Egyptian Government to pursue 
his Ph.D. studies at the University of Texas at Austin. 



References 

623 Baker, W.G., and Martyn, D.F. (1953). Electric Current in the Ionosphere, I, The 

624 Conductivity, Philos. Trans. R. Soc. London, Ser. A, 246, 281. 

625 Balsley, B. B. (1969), Some characteristics of non-twostream irregularities in the 

626 equatorial electrojet, J. Geophys. Res., 74(9), 2333-2347. 

627 Balsley, B. B., and D. T. Farley (1971), Radar studies of the equatorial electrojet at 
62a three frequencies, J. Geophys. Res., 76(34), 8341-8351. 

629 Balsley, B. B., and D. T. Farley, (1973), Radar studies of the equatorial electrojet at 

630 three frequencies, J. Geophys. Res. ,7 8, 7471. 

631 Benevolenskaya, E. E., and I.G. Kostuchenko, (2013), The Total Solar Irradiance, 

632 UV Emission and magnetic Flux during the Last Solar Cycle Minimum, Journal of 
Astrophysics, Vol. 2013, Article ID 368380. 

634 Bilitza D., and Reinisch, B. (2008), International Reference Ionosphere 2007: Improve- 

635 ments and new parameters, J. Adv. Space Res., 42, 4, pp. 599-609. 

636 Bowles, K.L., Balsley. B.B., Cohen, R. (1963). Field-aligned E region irregularities 

637 identified with acoustic plasma waves, J. Geophys. Res., 68, 2485. 

638 Buneman, O. (1963). Excitation of field-aligned sound waves by electron streams. Phys. 
Rev. Lett., 10, 285. 

640 Cohen, R., and Bowles, K.L. (1967). Secondary irregularities in the equatorial electro- 

641 jet. J. Geophys. Res. 72, 885. 



DRAFT 



January 14, 2015, 12:23pm 



DRAFT 



X - 34 EHAB HASSAN ET AL.: EQUATORIAL ELECTROJET TURBULENCE: BASELINE2D 

642 Farley, D.T. (1963). A plasma instability resulting in field-aligned irregularities in the 

643 ionosphere. J. Geophys. Res., 68, 6083. 

6M Farley, D. T., and B. B. Balsley (1973), Instabilities in the equatorial electrojet, J. 
Geophys. Res., 78,227. 

646 Farley, D.T. (2009). The equatorial E-region and its plasma instabilities: a tutorial. 
Ann. Geophys., 27, pp. 1509-1520. 

648 Fejer, B.G., Farley, D.T., Balsley, B.B., and Woodman, R.F. (1975). Vertical structure 

649 of the VHF backscattering region in the equatorial electrojet and the gradient drift 

650 instability. J. Geophys. Res., 80, 1313. 

651 Fejer, B. G., D. T. Farley, B. B. Balsley, and R. F. Woodman, (1976), Radar observa- 

652 tions of two-dimensional turbulence in the equatorial electrojet, J. Geophys. Res., 81, 

653 130. 

654 Fejer, B.G., and Kelley, M.C. (1980), "Ionospheric irregularities", Rev. of Geoph.and 
Space Phy, Vol.18, No.2, pp. 401-454. 

656 Forbes J.M. (1981). The Equatorial Electrojet, Rev. Geophy. and Space Phy, Vol 19, 
No. 3, pp. 469-504. 

ess Hanuise, C, and M. Crochet (1981), 5-50 m wavelength plasma instabihties in the 

659 equatorial electrojet, 2, Two-stream conditions, J. Geophys. Res., 86, 3567. 

660 Hoh, F.C. (1963). Instability of Penning-type discharge. Phys. Fluids, 6, 1184. 

661 Hysell, D.L., Drexler, J., Shume, E.B., Chau, J.L., Scipion, D.E., Vlasov, M., Cuevas, 

662 R., and Heinselman, C. (2007). Combined radar observations of equatorial electrojet 

663 irregularities at Jicamarca. Ann. Geophys., 25, 457473. 

664 Horton, W., and Hasegawa A. (1994), Quasi-two-dimensional dynamics of plasmas and 



DRAFT 



January 14, 2015, 12:23pm 



DRAFT 



EHAB HASSAN ET AL.: EQUATORIAL ELECTROJET TURBULENCE: BASELINE2D X - 35 

665 fluids, Chaos, pp. 227-251. 

eee Horton, W., and Ichikawa Y. H. (1996), Chaos and Structures in Nonlinear Plasmas, 

ee? World Scientific, Singapore. 

ees Kamenetskaya, G. Kh. (1971), Quasi-linear theory of the formation of inhomogeneities 

669 m the equatorial electrojet, Geomagn. Aeron., 11, 71-75. 

670 Kaw, P. K., (1972), Wave propagation effects on observation of irregularities in the 

671 equatorial electrojet, J. Geophys. Res., 77, 1323-1326. 

672 Knox, F. B., A contribution to the theory of the production of field-aligned ionization 

673 irregularities in the equatorial electrojet, J. Atmos. Terr. Phys., 26, 239-249, 1964. 

674 Kelley, M.C. (2009). The Earth's Ionosphere, Plasma Physics and Electrodynamics, 

675 2nd edition, Amsterdam ; Boston : Academic Press, 2009. 

676 Kudeki, E., D. T. Farley, and B. G. Fejer (1982), Long-wavelength irregularities in the 
e?? equatorial electrojet, Geophys. Res. Lett., 9, 684. 

678 Kudeki, E., D. T. Farley, and B. G. Fejer (1985), Theory of spectral asymmetries and 

679 nonlinear currents in the equatorial electrojet, J. Geophys. Res., 90,429. 

Kudeki, E., Fejer, B.G., Farley, D.T., and Hanuise, C. (1987). The CONDOR Equato- 

681 rial Electrojet Campaign: Radar results, J. Geophys. Res., 92, 13,561. 

682 Kudeki, E., and Farley, D.T. (1989). Aspect Sensitivity of Equatorial Electrojet 

683 Irregularities and Theoretical Implications, J. Geophys. Res., 94, Al, 426-434. 

684 Lu F., Farley, D.T., and W.E. Swartz (2008). Spread in aspect angles of equatorial E 
region irregularities, J. Geophys. Res., 113, A11309, doi:10.1029/2008JA013018. 

686 Machida, S., and C. K. Goertz (1988), Computer simulation of the Farley-Buneman 

687 instability and anomalous electron heating in the auroral ionosphere, J. Geophys. Res., 



DRAFT 



January 14, 2015, 12:23pm 



DRAFT 



X - 36 EHAB HASSAN ET AL.: EQUATORIAL ELECTROJET TURBULENCE: BASELINE2D 

93, 9993. 

ess Maeda, K., T. Tsuda, and H. Maeda, Theoretical interpretation of the equatorial 
ego sporadic E layers, Rep. Ionos. Space Res. Jap., 17, 3, 1963. 

691 Newman, A. L., and E. Ott (1981), Nonlinear simulations of type I irregularities in the 

692 equatorial electrojet, J. Geophys. Res., 86, 6879. 

693 Oppenheim, M. M., N. F. Otani, and C. Ronchi (1996), Saturation of the Farley- 

694 Buneman instability via nonlinear electron E x B drifts, J. Geophys. Res., 101. 

695 Oppenheim, M. M., and N. F. Otani (1996), Spectral characteristics of the Farley- 

696 Buneman instability: simulations versus observations, J. Geophys. Res., 101. 

eg? Oppenheim, M. M., Y. Dimant, and L. P. Dyrud (2008), Large-scale simulations of 

698 2-D fully kinetic Farely-Buneman turbulence, Ann. Geophys., 26, 543-553. 

699 Oppenheim, M. M. and Y. S. Dimant (2013), Kinetic simulations of 3-D Farley- 

700 Buneman turbulence and anomalous electron heating, J. Geophys. Res. (Space 
Physics), 118, 1306-1318, doi:10.1002/jgra.50196. 

702 Otani, N.F., Oppenheim, M.M., (1998). A saturation mechanism for the FarleyBune- 

703 man instability. Geophysical Research Letters 25, 18331836. 

70 4 Otani, N. F. and Oppenheim, M. (2006), Saturation of the Farley-Buneman in- 

705 stability via three-mode coupling, J. Geophys. Res. (Space Physics), 111, 3302, 
doi:10.1029/2005JA011215. 

707 Patra, A. K., Tiwari, D., Devasia, C. V., Pant, T. K., and Sridharan, R.: East-west 

708 asymmetries of the equatorial electrojet 8.3m type-2 echoes observed over Trivandrum 

709 and a possible explanation, J. Geophys. Res., 110, A11305, 2005. 

Pfaff, R. F., M. C. Kelley, E. Kudeki, B. G. Fejer, and K. Baker (1987a), Electric field 



DRAFT 



January 14, 2015, 12:23pm 



DRAFT 



EHAB HASSAN ET AL.: EQUATORIAL ELECTROJET TURBULENCE: BASELINE2D X - 37 

11 and plasma density measurements in the strongly driven daytime equatorial electrojet, 

12 1, The unstable layer and gradient drift waves, J. Geophys. Res., 92, 13578. 

Pfaff, R.F., Kelley, M.C., Kudeki, E., Fejer, B.G., and Baker, K.D. (1987b). Electric 
14 field and plasma density measurements in the strongly driven daytime equatorial 
is electrojet. 2 Two-stream waves, J. Geophys. Res., 92, 13,597. 

Pfaff Jr., R.F., Acua, M.H., Marionni, P.A., Trivedi, N.B., (1997a). DC polarization 
17 electric field, current density, and plasma density measurements in the daytime 
is equatorial electrojet. Geophysical Research Letters 24, 16671670. 

19 Pfaff Jr., R.F., Marionni, P. A., Swartz, W.E., (1997b). Wavevector observations of 

20 the two-stream instability in the daytime equatorial electrojet. Geophysical Research 
Letters 24, 16711674. 

22 Reid, G. C. (1968), The formation of small-scale irregularities in the ionosphere, J. 

Geophys. Res., 73(15), 1627-1646. 
24 Rogister, A., and N. D'Angelo (1970), Type II irregularities in the equatorial electrojet, 

J. Geophys. Res., 75(19), 3879-3887. 
26 Rogister, A. (1971), Nonlinear theory of type I irregularities in the equatorial electrojet, 

J. Geophys. Res., 76(31), 7754-7760. 

28 Rogister, A. (1972), Nonlinear theory of cross-field instability with application to the 

29 equatorial electrojet, J. Geophys. Res., 77 (16), 2975-2981. 

so Rogister, A., and E. Jamin (1975), Two-dimensional nonlinear processes associated 
si with "Type I" irregularities in the equatorial electrojet, J. Geophys. Res., 80, 1820-1828. 

32 Richmond, A.D. (1973). Equatorial electrojet. I. Development of a model including 

33 winds and instabilities. J. Atmos. Terr. Phys., 35, 1083. 



DRAFT 



January 14, 2015, 12:23pm 



DRAFT 



X - 38 EHAB HASSAN ET AL.: EQUATORIAL ELECTROJET TURBULENCE: BASELINE2D 

734 Rishbeth, H. (1997). The Ionospheric E-layer and F-layer dynamos - a tutorial review. 
J. of Atmos. and Solar- Terr. Phys. Vol.59, No. 15. 

Ronchi, C, Similon, P.L., and Sudan, R.N. (1989). J. Geophys. Res., Vol.94, No.A2, 
p.p.1317-1326. 

738 Sato, T., and T. Tsuda (1967), Computer study on nonlinear cross-field instability, 
Phys. Fluids, 10, 1262. 

740 Simon, A. (1963). Instability of a partially ionized plasma in crossed electric and 

741 magnetic fields. Phys. Fluids, 6, 382. 

742 Schmidt, M. J., and S. P. Gary (1973), Density gradients and the Farley-Buneman 

743 instability, J. Geophys. Res., 78, 8261. 

744 Skadron, G., and J. Weinstock, (1969), Nonlinear stabilization of a two-stream plasma 

745 instability in the ionosphere, J. Geophys. Res., 74, 5113-5126. 

746 Sudan R.N., Akinrimisi, J., and Farley, D.T. (1973). Generation of small-scale irregu- 

747 larities in the equatorial electrojet, J. Geophy. Res., Vol.78, No.l, pp. 240-248. 

748 Sudan, R . N. (1983a), Nonlineart heory of type I irregularities in the equatorial 

749 electrojet, J. Geophys. Res., 88, 40. 

750 Sudan R.N. (1983b). Unified theory of type-I and type-II irregularities in the equatorial 
electrojet, J. Geophy. Res., Vol.88, N0.A6, pp. 4853-4860. 

752 Sudan R.N., Gruzinov A.V., Horton W., and Kukharkin N. (1997). Convective 

753 turbulence in weakly ionized plasma, Physics Reports, Vol.283, No.l, pp. 95-119. 

754 Weinstock, J., and A. Sleeper (1972), Nonlinear saturation of 'type F irregularities in 

755 the equatorial electrojet, J. Geophys. Res., 77(19), 3621-3624. 

756 Wernik, A.W., Secan, J. A., Fremouw, E.J. (2003). "Ionospheric Irregularities and 



DRAFT 



January 14, 2015, 12:23pm 



DRAFT 



EHAB HASSAN ET AL.: EQUATORIAL ELECTRO JET TURBULENCE: BASELINE2D X - 39 

Scintillation", Adv. Space Res., Vol.31, No.4, pp.971-981. 

Woods, T.N., S. Cranmer, T. Hoeksema, and J. Khol (2010). "Irradiance Variations 
during this Solar Cycle Minimum," in SOHO-23: Underswtanding a Peculiar Solar 
Minimum, Eds., vol. 428 of ASP Conference Series, p. 63. 



DRAFT 



January 14, 2015, 12:23pm 



DRAFT 



X - 40 



EHAB HASSAN ET AL.: EQUATORIAL ELECTROJET TURBULENCE: BASELINE2D 



Table 1. Ionosphere Background Parameters 



Quantity 


Symbol 


Unit 


Solar Maximum 


Solar Minimum 


Electron Density 


N e 


m 3 


1.52 x 10 11 


1.60 x 10 9 


Electron Temperature 


T e 


°K 


192.1 


191.1 


Ion Temperature 


Ti 


°K 


192.1 


191.1 


Magnetic Field 


B 0 


T 


3.7 x 1(T 5 


3.5 x 10" 5 


Ion Acoustic Speed 


c s 


m/s 


318.5 


323.5 


Electron Thermal Speed 




m/s 


5.4 x 10 4 


5.4 x 10 4 


Ion Thermal Speed 


v u 


m/s 


225.2 


228.5 


Electron-Neutral Collision Frequency 




s- 1 


2.3 x 10 4 


2.4 x 10 4 


Electron-Neutral Collision Frequency 


^in 


s- 1 


3.2 x 10 3 


2.8 x 10 3 


Electron Gyro-Frequency 


Uce 


s- 1 


6.5 x 10 6 


6.1 x 10 6 


Ion Gyo-Frequency 


Uci 


s- 1 


113 


110 


Electron Larmor Radius 


Pe 


m 


8.3 x 1(T 3 


8.8 x 10" 3 


Ion Larmor Radius 


Pi 


m 


2 


2 
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Figure 1. The dynamic system geometry. The y-axis is parallel to the zonal direction and the z-axis is points 
upward. The magnetic field is out of the paper and parallel to the x-axis. The possible unstable linear-modes coupling 
with forwrad (a) or inverse (b) cascade energy transfer. 
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Figure 2. The growth-rate profile of the unified model as a function of altitude and horizontal wavenumber (k y ). 
It shows the dominance of Type-II instability at low altitudes up to the coupling region between Type-I and Type-II 
instabilities between 100-110 km. Above 110 km, Type-I instability dominates and Type-II instability disappears as a 
result of the reverse of the electron density scale length sign to be negative at that region. 
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Figure 3. A comaprison between the growth-rate, 7fc(fc H ), in the standard two-stream model (dashed-green) and 
the unified model with including both the ions viscosity and the electrons inertia (solid-black) , absence of ions viscosity 
(dashed-blue), and absence of electrons inertia (dashed-red) . 
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Figure 4. A comaprison between the real- frequency, u>k(k y ), profile of the standard two-stream model (dashed-green) 
and the unified model with including both the ions viscosity and the electrons inertia (solid-black). The unified model 
gives smaller phase velocity than the standard two-stream model. 



DRAFT 



January 14, 2015, 12:23pm 



DRAFT 



EHAB HASSAN ET AL.: EQUATORIAL ELECTRO JET TURBULENCE: BASELINE2D 



X - 45 



Spectrum of Density Perturbation 
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Figure 5. The Energy Spectrum of the Pertubed Electron Density at different states of the simualtion. The top 
panel shows the energy spectrum of the available modes accumulates at two peaks for two types of instabilities at low- 
and high-k values. The second panel shows the growth of modes at higher wavenumbers. The modes of k-values with 
negative growth-rate decay linearly. The third panel shows the state of energy distribution during the transition region. 
The fourth and the fifth panels show the redistribution of energy over the modes at all wavenumbers which verifies the 
nonlinear coupling between different modes in the system. 
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Figure 6. The Evoluation of the Pertubed Electron Density at t = 500 ms during the dominance of the linear 
terms on the simulation under the solar-minimum conditions. It shows the downward drifts of the high-density plasma 
and upward drifts of low-density plasma. 
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Figure 7. The Pertubed Electron Density at t = 600 ms during the transition from the linear to the saturation phase 
of the simulation under the solar-minimum conditions. It shows large size structures of plasma densities that are drifting 
horizontally and vertically at differnt velocities with maximum perturbation of 12%. 
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Figure 8. The Pertubed Electron Density at t = 1000 ms the saturation state of the simulation reaches 8% realtive 
to the background density under the solar-minimum conditions. It shows the very small structures embedded in the larger 
structures that are seen in the transition phase as a result of the coupling between Type-I and Type-II instabilities in the 
saturation state. 
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Figure 9. A close view of the large-scale structures formed during the transition phase (left-panel) 
and small-scale structures formed during the saturation phase of the simulation under the solar-maximum 
conditions. 
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Figure 10. A comparison between the the perturbed density maxima (top panel) and its stadard 
deviations (bottom panel) in the solar-maximum (blue-line) and solar-minimum (red-line) conditions. The 
standard deviations in the lower panel also defines clearly the growing, transion, and saturation phases 
during the simulation under two different solar conditions. 
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Figure 11. A comparison between the maximum perturbed electric field components in the horizontal 
(top panel) and vertical (bottom panel) directions in the solar-maximum (red-line) and solar-minimum 
(blue-line) conditions shows the effect of solar activity on the magnitude of the perturbed electric field 
components. 
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Figure 12. A comparison between the root-mean-square of the perturbed electric held components 
in the horizontal (top panel) and vertical (bottom panel) directions in the solar-maximum (blue-line) and 
solar-minimum (red-line) conditions. The dotted (dashed) lines show the root-mean-square of the positive 
(negative) component of the perturbed electric held. 
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Fi guxe 13. The inverse fourier transform of the particle density spectrum of all k y values at fc z =0 (m x ) for each 
time step during the saturation region. A phase velocity of 330 (m/s) can be estimated from the slope as the drift speed 
of the density to be close to the ion-acoustic speed at 105 m altitude under the solar-minimum conditions. 



DRAFT 



January 14, 2015, 12:23pm 



DRAFT 



